Mott Insulator to Superfluid transition in Bose-Bose mixtures in a two-dimensional 
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O 
(N 



(N 



M. Guglielmino/ V. Penna,^ and B. Capogrosso-Sansone"^ 

^ Dipartimento di Fisica, PoUtecnico di Torino, Corso Duca degli Abruzzi 24, 1-10129 Torino, Italy 

'^ Dipartimento di Fisica and CNISM Umta di Ricerca, 

PoUtecnico di Torino, Corso Duca degli Abruzzi 24, T10129 Torino, Italy 

'^Institute for Theoretical Atomic, Molecular and Optical Physics, 

Harvard-Smithsonian Center of Astrophysics, Cambridge, MA, 02138 

We perform a numeric study (Worm algorithm Monte Carlo simulations) of ultracold two- 
component bosons in two-dimensional optical lattices. We study how the Mott insulator to su- 
perfluid transition is affected by the presence of a second superfluid bosonic species. We find that, 
at fixed interspecies interaction, the upper and lower boundaries of the Mott lobe are differently 
modified. The lower boundary is strongly renormalized even for relatively low filling factor of the 
second component and moderate (interspecies) interaction. The upper boundary, instead, is affected 
only for large enough filling of the second component. Whereas boundaries are renormalized we find 
evidence of polaron-like excitations. Our results are of interest for current experimental setups. 
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In the last decade, a considerable amount of theoretical 
and experimental research has been devoted to the objec- 
tive of using ultracold lattice bosons and fermions to ad- 
dress many outstanding condensed matter problems via 
Hamiltonian modeling. Single species (Bose) Hubbard 
models, first introduced for fermions to describe electrons 
in solids and considered to be the minimal model for 
high temperature superconductivity, can be experimen- 
tally realized with ultracold atoms in optical lattices [1| 
and have been extensively studied 0, Q ■ If a second com- 
ponent is introduced, new fascinating phenomena and ex- 
otic quantum phases which cannot be accessed with sin- 
gle species atomic gases, become available. Prominent 
examples include the possibility of realizing quantum 
magnetic phases [J-Q ^ engineering ofFsitc interactions [a- 
[^ which would likely lead to supersolid states [3, [lO{ . 
studying disorder/impurities effects [ll|, [l2| • 

One of the fundamental questions to address in two- 
component systems is how the presence of a second 
species affects the Mott insulator (Ml)-superfluid (SF) 
transition. While MI phases of single component bosonic 
systems have been well established experimentally in 
one- [iJl , two- [lj| , and three-dimensions [l^ , the multi- 
component case results more challenging to understand 
as it involves additional degrees of freedom and hence a 
much richer physics. There have been a number of recent 
experiments on Bose-Bose piil - [l8l | and Fermi-Bose [2, Il9l - 
|21| mixtures which have investigated how the visibility of 
one component is affected by the presence of a second one 
upon varying relative densities of the two species and/or 
the interspecies scattering length. The overall finding is 
that the presence of a second component reduces the visi- 
bility of the first bosonic species. A number of theoretical 
explanations which include multi-band model [22|, self- 
trapping effect [12[ , changes in the chemical potential due 
to the presence of external harmonic confinement |23| . 
have been provided. 

In the present work we consider a homogeneous sys- 



tem of two-component bosons in a square lattice with 
repulsive interspecies interaction. We study how the MI- 
SF transition of the majority component is affected by 
the presence of a minority superfluid component upon 
varying the density of the latter and the interspecies in- 
teraction. This system can be realized by loading op- 
tical lattices with two different atomic species [la, Il7| . 
or the same atomic species in two different internal en- 
ergy states [l^ . Interspecies interaction strength Uab can 
be tuned either via Feshbach resonance or by changing 
the Wannier functions overlap (in the presence of state- 
dependent lattices). Intraspecies interactions Ua and Ub 
can also be tuned via Feshabch resonance. If the temper- 
ature is low enough, the system is accurately described 
by the two-component Bose-Hubbard Hamiltonian 
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-J J - "*j — rii'^j^]- Parame- 
ters ta, tb represent the hopping amplitudes, ^a, l^b the 
chemical potentials for the two bosonic species A and B, 
respectively. Index a = a,b identifies the two component. 
The phase diagram of model[T]is very reach. It includes 
several stable phases which vary from double MPs [2J], 
two independent SF's, and a mixture of one MI and 
one SF, to less trivialphases like supercountcrflow and 
checkerboard solid 0, [a] • In the following we present the 
first accurate results, based on path i nteg ral Monte Carlo 
simulations by the Worm algorithm [2^, for the MI-SF 
phase diagram of the majority component A. We study 
how the transition is affected by the presence of the mi- 
nority SF component B upon varying the density of the 



latter and the interspecies interaction. 

Our results can be summarized as follows. We find that 
the net effect of the SF minority component B is to screen 
the effective onsite repulsion of component A thus result- 
ing in a partial suppression of the Mott region, with the 
upper and lower boundaries of Mott lobes being very dif- 
ferently renormalized. Although the extension of the SF 
region would imply an increase of the system coherence 
and seems to contradict experimental findings, we point 
out that, in order to achieve a conclusive understanding, 
further investigation on, e. g., effects of the second com- 
ponent on the superfluid state of the majority species and 
of the presence of the harmonic confinement is needed. 
Finally, we find that, whereas boundaries are renormal- 
ized, particle and hole elementary excitations get dressed. 
The decay in imaginary time of the momentum-space 
Matsubara Green function of the component in the MI 
is reminiscent of the decay of the polaron Green function 
studied in the context of an electron coupled to optical 
phonons (see e. g. j26|). 

The large number of independent parameters occur- 
ring in Hamiltonian ([T]) involves an extremely structured 
phase diagram. In addition to the five parameters ta, 
Ua and Uab one must include the boson number of each 
species Na and Mb, controlled by chemical potentials 
fia and ^b- In the following we focus our attention 
on the first Mott lobe of species A, i.e. filling factor 
?1q = Afa/^I = 1, where M is the number of lattice sites. 
We use the MI regime of the single-component system, 
described by Hamiltonian ([1]) for Afb — 0, as the refer- 
ence case, and start our analysis by adding to the system 
a small amount of componen B {Afb << N'a) interact- 
ing with species A via a relatively weak interaction term 
Uab- We then proceed by progressively increase either 
Afb or Uab, while keeping fixed Ub = 10 if,, and tj, = ta 
(in the following we choose ta as energy unit). In order 
to obtain boundaries of the lobes, we either utilize the 
zero momentum Green function (GF) of component A 
(see [231 for details), or study na{fJ-a) curves at fixed Ub- 
The latter method is used whenever the imaginary time 
decay of the GF differs from what expected for quasipar- 
ticle/hole elementary excitation (see below). 

While in the absence of the second component the 
lobe boundaries are uniquely determined by varying pa- 
rameters ^a/Ua and ta/Ua, in the presence of the sec- 
ond species, the lobe boundaries are affected in different 
ways depending on the interplay of Uab with parameters 
of species A. We choose to use ta as our energy scale and 
determine the boundaries as Ua changes, while Uab is kept 
fixed for each lobe. We therefore expect the boundaries of 
the new lobe to coincide with those of the single-species 
case, i.e. Uab = or A/j, = 0, in the limit ta/Ua — >■ 0, since 
this limit also implies Uab/Ua -^ 0. We have performed 
simulations for linear system sizes L = 10, 20, 30, tem- 
perature T = ta/2L and have found no size effects within 
errorbars. The critical region at the tip of the lobes has 
not been studied as it is beyond the scope of this work. 

In Fig.[T]we show results referring to fixed Ub =0.1 and 
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FIG. 1: (Color online) First Mott lobe of species A with pa- 
rameters Ub = iota, tb = ta, rib = 0.1 and varying Uab = 10, 
15, 20ta, circles, squares, triangles respectively. The lobe for 
the single-species case, i.e. Afb — 0, is also plotted for refer- 
ence (diamonds), with data taken from [23|. Errorbars are 
within symbol size when not visible. Lines are to guide the 
eye. 



Uab = 10, 15, 2Qta corresponding to circles, squares and 
triangles, respectively. For comparison we also plot the 
single-species lobe (diamonds), taken from Ref [23|- As 
expected, all lobes tend to overlap to the single compo- 
nent one in the ta/Ua — >■ limit. In general, the change 
in shape of lobes involves a dramatic reduction of their 
extension even at low density of component B consid- 
ered here, Ub = 0.1. The suppression of the lobe be- 
comes more pronounced as the interspecies interaction is 
increased. The presence of species B therefore inhibits 
the formation of the species-A Mott state by favouring 
delocalization, resulting into a critical point ta/Ua renor- 
malized to lower values. 

Interestingly, while the lower boundary is strongly 
renormalized by the presence of component B, the upper 
boundary is essentially indistinguishable from the single- 
species case. Moreover, simulations done with Uab as 
large as ~ 40 ia suggest that, at nt ~ 0.1, phase separa- 
tion will occur before a sizeable shift in the upper bound- 
ary is detected. Indeed, in order to see any effect on the 
upper boundary, one has to go to larger density of com- 
ponent B (see Fig. [2] below). A simple physical argument 
can be made by considering the case where component B 
is reduced to just one particle. The energy, counted from 
chemical potential fia, required to create a hole excita- 
tion is reduced with respect to the single-component case 
as the hole represents a chance for the system to lower 
its energy: the empty site can be occupied by particle 
B. As a consequence, the interaction with species A is 
fundamentally turned off. The same effect cannot take 
place when a particle A is added to the lattice. In this 
case, due to the large dilution of species B which allows 
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FIG. 2: (Color online) First Mott lobe of species A with 
parameters Ub = Wta, tb = ta, Uab = 10 ia and varying 
Ub — 0.1, 0.5, 0.75, circles, squares, triangles respectively. The 
lobe for the single-species case, i.e. A/b = 0, is also plotted 
for reference (diamonds), with data taken from [2^ Errorbars 
are within symbol size when not visible. Lines are to guide 
the eye. 



the added particle A to avoid particle B, the energy cost 
for such an excitation remains substantially identical to 
that required when n^ = 0. These considerations can be 
extended to the case of Afb < < Ma and supported by the 
subsequent simple calculation. Denoting by E{Afa,M'b) 
the ground state energy of Hamiltonian ([T]) for negligible 
ta and tf), one finds the same particle-excitation energy 
gap A+ = EiMa + I, Mb) - E{MaMb) ^ Ua- ^ia for 
both Afb = and Afb ^ 0. Conversely, in the case of 
hole-excitation, where A_ = E{Afa - l,Afb)- E(Afa,Afb), 
one finds A_ = ^a for Afb = 0, and A_ ^ fJ-a — Uab < fJ-a 
for TVb ^ 0. If we refer to Fig. [T] (and Fig. O, with 
model parameters rescaled by I/C/q, we can observe that, 
for A/; ^ 0, A^/Ua -^ {/\-/Ua)M,=o when Ua -^ oo, 
confirming that different lobes will overlap in the limit 
ta/Ua — > 0. We can conclude that, at low B densities, 
the formation of a hole excitation is manifestly favoured 
by the presence of species B, whereas the mechanism for 
creating a particle excitation is essentially unchanged. 

Let us now turn to the discussion of the case of fixed 
Uab = iota and varying n^. Fig. [2] shows the resulting 
Mott lobes for Ub = 0.1, 0.5 and 0.75, circles, squares 
and triangles, respectively. The single component case 
is also plotted for reference (diamonds). As Ub is in- 
creased the lobe size is further reduced, although, even 
at large Ub = 0.75, the effect is less pronounced than 
what observed at Ub = 0.1 and large Uab ~ 15 — 20 to. 
The displacement of upper boundaries, though, becomes 
unequivocal for Ub = 0.5 and 0.75. Compared to the 
single component case, the energy required to create a 
particle excitation, counted from chemical potential /ia, 
significantly increases due to a large enough density of 



FIG. 3: (Color online) Imaginary-time decay of the zero- 
momentum GF of component A for L = lOT = ta/ 20 Ua = 

100 ia, Ub = 10 fa, tb =ta Ub ^ 0.1, Uab = 10 ia , Ma = 10 ia 

(triangles, lower x axis) and /ia = 89 ia (square, upper x axis) 
for hole and particle excitations, respectively. Solid lines are a 
linear fit. The linear decay of In G(0, r) on the hole excitation 
side is recovered only at later imaginary times, implying Z- < 
1 (polaron-like behavior). On the particle side, instead, the 
expected linear behavior is observed. 



component B. The latter discourages double occupancy 
of species A due to interaction Uab ■ 

In this scenario, a discussion of the imaginary time de- 
cay of the momentum space GF is appropriate. The GF 
is defined as follows: 



G(z,r) = (r.at(r)ao(0)) 



(2) 



where r is the imaginary time, and 7^ represents the 
time ordering operator. Using Lehman expansion and 
extrapolation to the t — > ±oo limit one readily finds: 
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The two limits describe single-particlc/hole excitations 
in the MI phase. Here Z± and e± are the particle/hole 
spectral weight (or Z-factors) and energy, respectively. 
In the grand canonical ensemble, excitation energies are 
measured relative to the chemical potential. Chemi- 
cal potentials ^± for which the energy gap for creating 
the particle/hole excitation with p = vanishes, are 
then retrieved by the exponential decay of G{Q,t). If 
elementary excitations arc interpreted as polarons. the 
spectral weights Z± carry information on which frac- 
tion of the bare-(quasi)particle state is present in the 
polaron cigenstate. In the polaron GF, the exponential 
decay is expected to be recovered only at later imaginary 
times, therefore, (properly normalized) spectral weights 
are smaller than 1, Z± < 1 [2y|. This is precisely what 
we have observed. In Fig. [3] we plot Ga{0, t) in log scale, 

for Ua = lOOta, Ub = iota, tb = ta, Ub = 0.1, Uab = 10 ta. 



jia ~ 10 ta and fj,a = 89 ta corresponding to hole (trian- 
gles) and particle (squares) excitations respectively. By 
using the periodicity of GF in r, we show the plot for 
positive values r g [0,/3], /3 = 1/T. Since we also per- 
form the transformation r — > — t, the r interval relevant 
to hole (particle) excitations is that starting from t ~0 
(ending at t = /3). Solid lines are linear fits, with the 
y-intcrsect providing spectral weights Z±. Clearly, the 
decay on the hole excitation side initially differs from the 
purely exponential one expected. The latter is recov- 
ered at later imaginary times, implying Z_ < 1. This is 
the typical behavior of the polaron GF. On the particles 
side, instead, the expected exponential decay is observed 
(here the lobe boundary coincides with the single com- 
ponent case). We notice that, wherever lobe boundaries 
differ from the single species case, we always observe a 
polaron-likc behavior of the GF, both for particle and 
hole excitations. We interpret this observation as dress- 
ing of elementary excitations, which we ascribe to the 
coupling of quasiparticles and holes of species A with 
quasimomcntum modes of SF species B. The latter play 
the role of phonon modes within the tight-binding picture 
of Frohlich's model [2^. A simple argument seems to 
confirm this picture. Calling Rj the position of j-th site 
and expressing species-B operators in terms of momen- 
tum modes bj = ^^^b^ exp{i'k ■ Rj) / y/M we can think 
the zero-momentum mode to be macroscopically occu- 
pied in the SF regime, which allows us to make use of 
the Bogoliubov approximation 6o = ^J — VA/fc. Hamilto- 
nian ([ij reduces to H ^ Hu ~ flaNa — fJ,bNb + Hpoi, where 
the Hu collects all the intraspecies interaction terms, and 
i^a = l^a~ Uab '^b IS a reuormalized chemical potential for 
species A. Hpoi reads 
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where Wk = 2tb [cos{dkx) + cos{dky)] is the dispersion 
relation for species-B atoms, d is the lattice spacing, 
and Uab = UabVT^/M. The last term in (|j) is of 
particular interest since it reproduces the interaction 
scheme between phonons and electrons in the conductive 
band leading to polaron excitations [2^. We notice that 
polaron physics in the context of impurities embedded 
in ultracold atomic systems has been a topic of great 
interests in the last fewyears, both experimentally and 
theoretically (see e.g. 0, [l^, [2^, |30|)- In the present 
work, though, polaron physics emerge not from dressing 
of impurities but of elementary excitations. 

Concluding, we have performed a numeric study of the 
MI-SF transition in the presence of a second SF bosonic 
species showing that, at fixed interspecies interaction, 
the upper and lower boundaries of the Mott lobe are 
very differently modified. The overall effect of the 
second component is reducing the extension of the MI 
lobe. We have also found that, whereas boundaries are 
renormalized, the decay of the GF exhibits a polaron-like 
behavior. In the future it would be interesting to study 
how the second component affects the superfluidity of 
the majority component in the presence of the harmonic 
confinement, as directly relevant to experiments. The 
investigation of the effective mass and the dispersion 
relation of polaronic excitations at strong, weak and 
intermediate coupling regimes also deserves further 
study. These aspects will be considered elsewhere |3]| . 
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